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CALCULATION OF THE LATERAL STABILITY OF A 
DIRECTLY COUPLED TANDEM-TOWED FIGHTER AIRPLANE 
AND CORRELATION WITH EXPERIMENTAL DATA 
By Robert E. Shanks 


SUMMARY 


A theoretical method is presented for predicting the dynamic lateral 
stability characteristics of an airplane towed in tandem by a much larger 
airplane. Values of period and time to damp to one-half amplitude and 
rolling motions calculated by an analog computer have been correlated with 
results of two experimental investigations conducted in the Langley free- 
f light tunnel which were part of a U. S. Air Force program (Project FICCN) 
to develop a satisfactory arrangement by which a bomber could tow a para- 
site fighter. In general, the theoretical results agree with the exper- 
imental results. 


INTRODUCTION 


As part of a U. S. Air Force program (Project FICCN) to develop a 
satisfactory arrangement by which a bomber could tow a parasite fighter, 
two experimental investigations have been conducted in the Langley free- 
flight tunnel to determine the dynamic stability characteristics of models 
of two types of fighter airplanes towed by a bomber. On the basis of the 
results of these investigations reported in references 1 and 2 and of sub- 
sequent full-scale flight tests, a satisfactory tandem directly coupled 
configuration was selected as the final arrangement. In order to provide 
a means for predicting the lateral stability characteristics of future 
directly coupled aircraft configurations of this type, the equations of 
lateral motion are derived for the towed airplane. The period and damping 
have been calculated from the characteristic modes of the system and time 
histories of the rolling motions in response to specified aileron disturb- 
ances were also obtained by using this method for the conditions of the 
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model tests of the Republic F-84E airplane and the Chance Vought 
F7U-3 airplane of references 1 and 2, respectively. The results of the 
calculations are correlated with the experimental data. 


SYMBOLS 


All forces and moments are referred to the conventional stability 
system of axes (see fig. l) . 


W 

m 

S 

K 

b 

t 

V 

v 

<1 

P 

P 

n 


weight of model, lb 

mass of model, slugs 

wing area, sq. ft 

spring constant, ft-lb /radian 

wing span, ft 

time, sec 

airspeed, ft/sec 

sideslip velocity, ft/sec 

dynamic pressure, i p y , lb/sq. ft 

air density, slugs/cu ft 
relative-density factor, m/pSb 

angle of attack of principal longitudinal axis of inertia, 
deg (see fig. l) 

radius of gyration about principal longitudinal axis of 
inertia, ft 

radius of gyration about principal normal axis of inertia, ft 


/ ^ p 2 2 

radius of gyration about X-axis, rx Q 'cos tj + kg sin tj, ft 


cos t) + ky sin t), ft 
u •*’o 


radius of gyration about Z-axis, 
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*xz 



product-of-inertia factor. 




sin t| cos t] 


L b 


K XZ = 


^XZ 

2 


3 


x horizontal distance from, tow-attachment coupling to center 

of gravity measured along longitudinal stability axis; 
values of x must be determined for each angle of attack, 
ft (see fig. 1) 

z vertical distance from tow-attachment coupling to center of 

gravity measured along normal stability axis; values of z 
must be determined for each angle of attack, ft (see fig. l) 


A,B,C,D,E 

A 

D 

s 

P 

T l/2 


coefficients of the stability characteristic equation 
root of stability characteristic equation 
differential operator, d/ds 

nondimensional time parameter based on span, tV/b 
period of lateral oscillation, sec 
time to damp to one-half amplitude, sec 


0 angle of sideslip, radians 

i|r angle of yaw, radians 

(f angle of bank, radians 

r yawing angular velocity, di|f/dt, radians/sec 

p rolling angular velocity, d0/dt, radians/sec 

Y side force, lb 

N yawing moment, ft- lb 

AnFSfflBSSiS' 
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rolling moment, ft- lb 

BY 

SP 

Bn 

Bp 

Bl 

Sp 

Bl 

Bt 

Bn 

B^ 

Bl 

Bt 

By 

Bp 

Bn 

Bp 

BL 

Sp 


lateral- force coefficient, 
yawing-moment coefficient, 
rolling-moment coefficient. 


Y/qS 

N/qSb 

L/qSb 
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CL = 


Hp 


_ aC n 


d 2V 


Sc 

8 


l 

pb 

2V 


= 


SCy 

d 2V 


Jn r 


C, = 


52n 

d 2V 


SC 

s 


rb 

2V 


°yo = 


SCy 

Jb 

a 2? 


J n- 


*0 


SCn 

(3b 

S 2V 

5£i 

Jb 

a 5v 


c w 


weight coefficient, W/qS 


-K 


spring-constant coefficient, K/pSV b 


Subscripts : 
c constraint 


o 


control or initial 
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THEORETICAL METHOD 
Assumptions of the Analysis 


In addition to the customary assumptions of small disturbances and 
linear derivatives, the following assumptions were also made: 

The motions of the towing airplane arising from the fighter motions 
are assumed to be negligible because of the large difference in the 
relative size of the two airplanes, and the towing airplane is assumed 
to remain in steady level flight. 

The fighter airplane is assumed to be rigid, and it is assumed that 
the tow-attachment coupling is mounted to the towing airplane without 
flexibility; therefore, elastic effects may be ignored. 

The air flow in the vicinity of the fighter is assumed to be smooth 
and unaffected by the towing airplane. 


! a 


1 


Equations of Motion 

The equations of motion of the towed airplane were developed by 
using Lagrange's method of undetermined multipliers as a means of taking 
account of the constraint on the airplane. This method is based on 
d'Alembert's principle and is discussed in detail in reference 3. 

The equation of condition which defines the sidewise velocity at the 
constraint or tow-attachment coupling is 

v c = v + Vt + xt^ - = 0 (l) 

Solving for the sidewise velocity of the airplane center of gravity from 
equation (l) results in 

v = z0 - xi(f - (2) 

and 

V = 20 - jojf - V+ (3) 

The lateral equations of the free airplane are: 

Side-force equation, 

m(v + Vi) = Yp | + Y p | + W0 + Y Q = Y (k) 
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Yawing-moment equation, 

V - W “ M |j | + V + y + 2 + M 0 = H (5) 

Rolling-moment equation, 

Z X0 - X XZ^ = L P ^ + kty 1 ^ + L 0$ + L P ^ + K00 + L 0 = L (6) 

Use of d'Alembert's principle and the Lagrange multipliers res ul ts in 
the following equation: 


j^m(v + Vty) - Y~|sy + [l z t - " N ]sty + 

j^ljr0 - lyy.ty - lJ 50 + A|j5y + x Sty - z S0j = 0 (7 a ) 

or 

£m(v + Vty) - Y + A]5y + |jE z ty - Ixz0 - N + xA^Sty + 

j^I x 0 - “ L " zAj&0 = 0 (7b) 

Because of the inclusion of A, the virtual displacements 5y, Sty, 
and 50 are arbitrary and, therefore, the coefficients of each displace- 
ment must vanish. Then, solving for A from the coefficient of 5y 
yields 

A = Y - m(v + Vty) (8) 

Substituting for A and for v and v from equations (2) and (3) in 
the yawing-moment and rolling-moment equations gives 


I z ty - Ixz0 - f N p( z 0 - - v ^) - W ty^ - - i K p ( Z 0 - xty - Vty) + 

i Y p (z0 - xty - Vty) + i Y p (z0 - xty - Vty) + Y^ty + Y00 + ¥0 - 


m(z0 - xty - Vty) - mVty 


= N 0 - xY 0 


(9) 




hff 

I 


m 
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I X 0 - IxZ^‘ " ^ L (i ( z 0 " 4 - Vt) - Lji - - i L p (z0 - xt - Vt) - 

K00 - z(i Y£(z0 - x¥ - W) + i Y p (z0 - xi|f - Vi|f) + + Y00 + 

W0 - m(z0 - xi[r - Vijr) - mVfJ = L 0 + zY Q (10) 

These equations may he expressed in nondimensional form by dividing each 

I O 

by ipSV b and by using the following relationships and the operator D 


m 

pSb 


Ky = — 


D = A 
ds 


v _ ^Z 

KZ - “ 


ds = — dt 
b 


K 


XZ 


xz 


The equations in nondimensional form are: 
Yawing-moment equation, 

2 . ,^2 n v _ iM 2 , lj2_. 


2hK 7 + &(*) + i £ Cn. _ 14 s ) (V. D> + i C n . - 
_ L \b) 2b p 2 lb/ ip (_2 n p 


1 — Cy. - i C n + 

2 b x p 2 r 


1 x Cy + £ C n - 

2 b b P Vb 


+ ( Cn p “ b ^p)^ “ ( 2(iKxz + 2|a 


xz 


1 z c 1 xz Gy^D 2 0 - C n -iS(L + z C n - 2“ (V \D0 

2 b n p 2 b 2 ^pj ^ \2 n P 2 b i p b n p b 2 ^p/ y 


l C W)0 = C no - * ^ 


( 11 ) 
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Rolling-moment equation. 



2 i iSZ + lxc 7 .+i^.G v \D 2 t + fl C,. + i S. Gy. 
b 2 2 b 2 P 2^2 \2 Ip 



1 Z p 

2b% 






+ 



- ~ C7 • 

2 1 4 P 




1 z p 

2 b % 



+ 




( 12 ) 


When 0o e ^ S is substituted for 0 and \ir 0 e is substituted for »|r 

in the equations written in determinant form, A must be a root of the 
stability equation: 


k 3 2 

AA + BA + CA + DA + E = 0 


(13) 


where 


A = V 


2 2 /„\ 2 2 /.A 2 2 2 „„ 

K X % + (§) Kz + \p) K X " K XZ ~ 2 — K XZ 




b " (bf K x 2 % “ f " (|) *Z % + 


,2 2 


X t/- rt Z. v p _i_o 2CZ v n 

b K XZ Cz (3 " b K XZ Cn (3 2 -2 
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B = 2|i 


- (f) 2 ^ - i *xz Cn e + f %%, - ffiVor, + 


£ K^C, + 2^2. C v 
t 'XZ l o 2 X 


+ H 


2, /„\2 

'*p vw ’t? 


- g % 


§ + f - KxV - ft) S- - Wi r - rl c !r + f *x 2 °r r 


5 K xz% + *xS + (§f% + Kxz ° ! p + P % - b *x\ + f *bi% : 


i - 2 C n .C 7 + 2 C 7 .Cn + Cy.C 7 - fe') Cy.Cy +#07.0!, 
4|_ b n p Z P b z p “p W ^p Z P W z p^p 13 l p “r 


t 0 7 C n . + Cy.C n - O.Cy- - fi.Cy - 2% C 7 Cy 

+ z ~ n p W Y P n r \b/ n p^r b 2 n p^p b 2 


-Cy.Cn +2 f C C z - 
b 2 P P b 2 P r 


C _ 2p -K z 2 ^ C W + C K j + | K XZ C W - ^ C K + % 2 C n ^ + C n ^ + K xz C 2 ^ + 

* * ■Witp] + I [- f %(f °w + c k) - ( 5)S°» + 


^ C z - S K 2 Cy + £ 

b 2 ‘-p b A x n >, 

,2 


xz 


^ S c x 


b “p *p 
\2 


#) C C K - S fi c 7 + f£) C n + £ C, C n + ^ C. C. 

b/ X R K -h n R i-n lb/ 1 b z„ T) . P x n ■ 


P TP 


‘•p “P b 2 "P “P 


P Y P W 'p Y P b ^p^ r ^p n r b ~ n p 


c z„ c y_ + r C z C. + # Cy L - # C n C. + 


S| Cy C - *% C % - 

b 2 P l r b 2 p^r Vb 


z ' 2 c. 


’ n Q Cy 

P x r 


+ i 


— C„ C 7 + — Cy Cj + 

n p z p b ^p Zp 


C n 07 - — Cy Ci + C 7 .C,, — C 7 Cj 7 - — Cy Cj! + — Cy . Cri - 

r *p b x r 4 p 4 p Hp l T °p b P b x p n p 


P P 


b S^P + b Cn r° Y p ” b Cz 6 Cy p + t C Z r C Y ; 


P A P 


t r -*-p 
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D = 







P 6 P 


c n r, c l„ + £ C Y r C Z t 


L P "p 


+ C-, C 


n 


P 




P s 


- £ c. 

b 


n P% 


- c 7 

b 1 


Cy 

P ^P 


E = “ Cn p(t °w + C k) - f ( C Z p C W - °r p C K ) 


The following relationships were used to determine the period and time 
to damp to one-half amplitude from the roots of equation (13) : 


P = 2jtb 
cnV 


T l/2 = 


0.693b 

aV 


where a and o> are the real and imaginary parts of the complex 
roots A = a t io}. 


Calculations of Period and Damping 


The period and time to damp to one-half amplitude of the lateral 
oscillations for the models of the F-84E and F7U-3 airplanes were 
calculated for the conditions listed in table I. The conventional 
static-lateral-stability derivatives used in equations (ll) and (12) 
were determined from force tests in the Langley free-flight tunnel and 
the conventional rotary lateral derivatives used in these equations were 
estimated by methods presented in reference 4. The values for the 


moment arms 
lowing table: 




used in these equations are given in the fol- 
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t 

& 


• ••• 

• • • 
• • • 

> • •• 

* • • 

J* *• 

■ •••• 

• • 
* 

&•••• 

••• 

»■ 


Airplane 

Loading 

x/b 

z/b 


I 

0.399 

- 0.076 

F-84E 

II 

.398 

-.093 


III 

.397 

-.111 

F7U-3 


.488 

-.151 


For both models, calculations were made for systematic changes in spring 
restraint in roll corresponding to the springs used in the experimental 
investigations. For simplicity in discussing the results, the springs 
were numbered from 1 to 5 and the following table gives the values of 
the various spring parameters: 

>1 


Spring 

Full-scale 
spring constant, 
ft-lb/radian 

Value of 

C K + § °w 


F-84E loading 

F7U-3 

I 

II 

ill 

None 

0 

- 0.0388 

-0.049 

- 0.058 

-0.0917 

1 

100,000 

-.1358 

-.146 

-.155 


2 

180,000 




- .1762 

3 

200,000 

-.2448 

-.255 

-.264 


4 

235,000 

| 



-.2019 

5 

490,000 




-.3248 


The spring-restraint parameter 


C K + § C W 


depends not only on the 


spring constant but also on the weight of the airplane. In the pre- 
ceding table, therefore, for the F-84E model there are three values of 

n + 'L a T for each spring corresponding to the weight at each loading 
K b w 


condition. 
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For the F-84-E model, calculations of the period and time to damp 
to one-half amplitude were also made to show the effect of artificial 
increases in the damping in yaw for loadings I, II, and III (p = 18 . 36 , 
23*0, and 27 . 3 > respectively) with no roll spring restraint and for the 
loading III, spring 3 case. 

The values of the artificial damping in yaw for the model of the 
F-84E airplane were calculated from the measured response characteristics 
of the rate-gyro — servomechanism system and the measured rudder effec- 
tiveness. The yaw-damper conditions of the tests were represented as 
nearly as possible in the calculations by the three different values of 
artificial damping in yaw which are listed for the model of the F-84E in 
table I. It should be noted that these values are considerably larger 
than the first estimates reported in reference 1, which were in error. 


Calculation of the Motions 

Motions were calculated with the Reeves Electronic Analog Computer 
(REAC) by using equations (ll) and (12) for all the conditions summarized 
in table I. The disturbances used for these calculations represented the 
disturbances measured from the test records as nearly as possible. The 
duration and sequence of the aileron controls used for the calculations 
corresponded to those measured from the test records. The aileron rolling 
moments were obtained from force tests of the model. 

In the tests of reference 1 for the stable conditions, the 
F-84E model was disturbed by an abrupt 35° aileron deflection and then 
the controls were centered for the rest of the test. This disturbance 
was represented for the REAC calculations by applying a rolling moment 
corresponding to 35° aileron deflection (^l Q = 0 * 05 j for 0*2 second. 

The time of 0.2 second was an approximate value obtained by averaging 
the values read from the film records for several cases. For the 
unstable test conditions, the oscillations were started by a gust in the 
tunnel airstream after the model had been steadied with the controls or 
with a launching line and then released with the controls fixed for the 
test record. For the calculations involving unstable oscillations, the 
gust was represented by a rolling moment one-tenth as large as that used 
for the stable cases = 0 . 005 ^ . 

In the tests of reference 2 for the stable configurations, the oscil- 
lations of the F7U-3 were started by alternate left and right aileron 
controls in phase with the rolling displacement to produce large- 
amplitude motions. Because the control sequences differed somewhat for 
the various cases, the control direction and duration were obtained from 
the film records for each condition and duplicated for the REAC calcula- 
tions. Rolling and yawing moments corresponding to 30° aileron deflection 
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fCi Q = 0.024, C no + ^ Cy o = 0.002^ were used for these cases. For the 
unstable cases, a single rolling moment = 0.02^ was applied for 


0.1 second to represent the tunnel gust. 


A preliminary comparison of the experimental data with the motions 
calculated as described in the preceding paragraphs showed that in most 
cases the calculated initial response was different from that obtained 
in the tests, apparently because the disturbance used in the calcula- 
tions did not exactly correspond to that used in the tests. The ampli- 
tudes of all the calculated motions were therefore adjusted so that the 
curve at the first peak following the disturbance (called time 0) 
matched exactly the experimental records. Mathematically this was 
equivalent to adjusting the magnitude, but not the duration, of the 
disturbance. It was felt that this procedure was justified because the 
primary interest in this investigation was not the initial response to 
the disturbances but rather the motions following the initial response. 


RESULTS AND DISCUSSION 


The results of the calculations and the corresponding data obtained 
from the experiments are presented in figures 2 to 9* The calculated and 
experimental values of period and damping are compared in figures 2 and 7 
for systematic variations of roll spring restraint ^as measured by 

Cr + Cy^) for the F-84E and F7U-3 models, respectively. The experi- 
mental values of time to damp to one-half amplitude were obtained from 
test records of references 1 and 2 by taking the average values from the 
damping envelopes of the first several cycles. The calculated and exper- 
imental rolling motions for the same conditions are compared in figure 3 
for the F-84E model and in figures 8 and 9 f or the F7U-3 model. In the 
case of all the motions calculated on the REAC, only the rolling motions 
are presented because only the rolling motions were obtained from the 
tests of references 1 and 2. 

The effect of the yaw damper for various loadings of the F84E model 
is shown in figure 4 by both the experimental and calculated data. The 
calculated and experimental motions are compared for the same conditions 
in figures 5 and 6. 


F-84E Results 

Effect of roll spring .- For all the F-84E conditions without the yaw 
damper, the calculations of period and time to damp to one-half amplitude 
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predict two oscillations which will be designated as the more stable 
oscillation and the less stable oscillation. Comparison of the calcu- 
lated and experimental damping values in figure 2 indicates that the 
oscillation obtained in the tests corresponds to the less stable mode 
predicted by the calculations, although in a few cases the experimental 
periods appear to agree more closely with the periods calculated for the 
more stable oscillation. In general the calculated values of period 
and time to damp agree with the experimental values for all the condi- 
tions shown in figure 2. No experimental value for damping is presented 
for the loading I, roll spring 3 case because the irregular motion made 
it impossible to determine the time to damp. The comparison of the test 
and calculated motions in figure 3> however, shows that the calculations 
do predict the type of motion for this condition. The slight irregu- 
larity in the test record of the loading I, roll spring 1 case in fig- 
ure 3 appears to be evidence of a short-period oscillation superimposed 
on the longer period motion but it is not apparent in the calculated 
motion. 

Figure 2 shows that the calculated periods of the two oscillations 
are usually about the same, except for the loading I case where they 
become increasingly different with increasing roll spring restraint. 

This large difference in the periods and the larger initial response 
of the more stable mode to the disturbance may explain the irregular 
motion of the loading I, roll spring 3 case. The two modes get out of 
phase early in the record, and when the amplitude of the stable mode is 
about the same as that of the less stable mode, the combination produces 
an unevenly damped motion before the more stable mode disappears. This 
irregular motion cannot be attributed to the type of control disturbance 
used because a smooth motion was obtained in the loading II, roll 
spring 3 case with the same disturbance. 

Effect of yaw damper .- The calculations for period and time to damp 
for the F-84E model with the yaw damper operating predict one oscilla- 
tion and two very stable aperiodic modes. Figure 4 shows the effect of 
the yaw damper on the period and damping of the oscillation for the 
three loading conditions corresponding to values of p of 18.36 
(loading I), 23.0 (loading II), and 27.3 (loading III). Only the values 
of period and time to damp for the predominant, less stable oscillations 
of the cases without the yaw damper are presented in figure 4 for com- 
parison with the values for the yaw-damper cases. The calculated aperi- 
odic damping values for the yaw-damper cases are presented in table I. 
The calculated and experimental motions showing the effect of the yaw 
damper are compared in figure 5 Tor the three loading conditions with- 
out roll spring restraint and in figure 6 for the loading III, roll 
spring 3 condition. 

The calculations predict the effect of the yaw damper fairly well, 
although in some cases the quantitative agreement between the test data 
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and the calculations is not so good. The test and calculated periods 
shown in figure 4 are only in fair quantitative agreement for loadings II 
and III but are in qualitative agreement in that both increased with 
increasing relative-density factor. The test and the calculated damping 
results are also in qualitative agreement in showing that increases in 
damping are provided by the yaw damper. Most of the quantitative differ- 
ence results from the differences between the calculated results and the 
test points for the no-yaw-damper cases. The calculated and experimental 
motions in figures 5 and 6 are in qualitative agreement, and the chief 
differences are in the periods of the calculated and experimental records. 


F7U-3 Results 

All the calculations made for the FJU-3 model predict two oscilla- 
tions. The calculated and experimental values for period and time to 
damp are compared in figure 7» Experimental values for the period are 
presented for all the conditions, but experimental damping values are 
presented only for the no-spring condition and for the roll spring 2 
case because the irregularity of the motions for the spring 4 and 5 
cases makes it very difficult to determine the damping values . The 
results in figure 7 show that the model had somewhat better damping 
than that predicted by the calculations for the no-spring condition 
and for the roll spring 2 case. 

The calculated and experimental records for the four spring- 
restraint conditions are presented in figure 8. As indicated by the 
data of figure 7j "the experimental motions show more damping than the 
calculated motions for the no-spring condition and for the roll spring 2 
case. The comparison of the records for the roll spring 4 and 5 cases 
shows fairly good agreement between the experimental and calculated 
motions. In the roll spring 4 case, the calculations predict the same 
irregular motion in the first three cycles that was found in the exper- 
imental records. 

The data of figure 7 show that the periods of the two oscillations 
are close for the low-spring-restraint conditions but the periods become 
increasingly different with increasing spring restraint and, as in the 
case of the loading I condition for the F-84E model, both the calculated 
and experimental motions of the F7U-3 model in figure 8 show uneven 
damping in the first several cycles with roll springs 4 and 5 where the 
periods of the two predicted modes were most different. 

The disturbance sequences were somewhat different for each of the 
F7U-3 model tests shown in figure 8. In order to verify the results of 
the F-84E tests which indicated that the type of disturbance was not 
responsible for the irregular motion, calculations were made for the 
roll spring 2 case by using the same disturbance as that for the roll 
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spring 4 case. Results of these calculations are presented in figure 9> 
together with results for the roll spring 4 case. The data of fig- 
ure 9 show that the oscillation for the roll spring 2 case is still 
smoothly damped with this disturbance. It appears, therefore, that, as 
in the case of the F-84E, the type of disturbance is not the cause of 
the irregular damping. 

The residual motions in figure 8 from about 5 seconds on for the 
roll spring 4 case and from about 2.5 seconds on for the roll spring 5 
case are evidently not oscillations predicted by the calculations. 

Since the calculations predict good damping for all the modes of these 
cases, it appears that the residual motions are probably caused by tur- 
bulence in the tunnel airstream as suggested in reference 2. 


CONCLUSIONS 


A theoretical method is presented for predicting the dynamic lateral 
stability characteristics of an airplane towed in tandem by a much larger 
airplane. Calculated values of period and time to damp to one-half ampli' 
tude and calculated rolling motions have been correlated with results of 
two experimental investigations conducted in the Langley free-flight 
tunnel. In general, the theoretical results agree with the experimental 
results. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 31> 1955* 


Approved : 


Robert E. Shanks 
Aeronautical Research Scientist 


dh 
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